Entanglement and sensitivity in precision measurements 
with states of a fluctuating number of particles 



P. Hyllus 1 , L. Pezze 2 and A. Smerzi 1 
1 INO-CNR BEC Center and Dipartimento di Fisica, Universitd di Trento, 1-38123 Povo, Italy 
2 Laboratoire Charles Fabry de Unstitut dOptique, CNRS and Univ. Paris-Sud, 
Campus Poly technique, RD 128, F-91127 Palaiseau cedex, France 
(Dated: October 13, 2010) 

The concepts of separability, entanglement, spin-squeezing and Heisenberg limit are central in 
the theory of quantum enhanced metrology. In the current literature, these are well established 
only in the case of linear interferometers operating with input quantum states of a known fixed 
number of particles. This manuscript generalizes these concepts and extends the quantum phase 
estimation theory by taking into account classical and quantum fluctuations of the particle number. 
Our analysis concerns most of the current experiments on precision measurements where the number 
of particles is known only in average. 
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Introduction. After the first pioneer works of the 
80's and early 90's [H-Q, the field of quantum-enhanced 
metrology is witnessing a revival of interest due to pos- 
sible implications on fundamental questions of quantum 
information and technological applications. Recent theo- 
retical analyses are mainly focusing on the interplay be- 
tween spin-squeezing, entanglement and phase sensitiv- 
ity 0, Q , the set of fundamental sensitivity bounds @- 
P and the quest for optimal phase estimation protocols 
[lol [IH . This activity has came along with breakthrough 
experimental investigations with photons ions [l3| . 
cold atoms [13] and Bose-Einstein condensates (BECs) 



can be written as 



15l |16| . The interplay between theory and experiments 
is playing a crucial role for the development of the field. 
It should be noticed, however, that while the theoret- 
ical investigations have been developed in the context 
of systems having a fixed, known, number of particles, 
most experiments have been performed with a fluctuating 
number of particles. The consequences of classical and 
quantum fluctuations have been generally disregarded. 



The goal of this manuscript is to extend along this 
direction the quantum metrology theory, and, in par- 
ticular, to discuss the relation between separability, en- 
tanglement, spin- squeezing and sub shot-noise sensitivity, 
and to settle the concept of the Heisenberg limit. In this 
framework, we recognize as entangled, for instance, the 
state suggested by Caves in '81 in the context of Mach- 
Zchnder interferometry [l|, which provided the first pro- 
posal for sub shot-noise phase estimation. 

Separability and Entanglement. The fluctuation of the 
total number of particles can have a classical and/or 
quantum nature. It is generally believed that quantum 
coherences between states of different numbers of parti- 
cles do not play any observable role because of the exis- 
tence of superselection rules (SSR) for the total number 
of particles 17 -III. Therefore, the most general states 
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where p( N > is a state of N particles. We extend the usual 
definition of separability of states of N particles [2(| by 
considering such incoherent mixtures as separable if they 
can be written as 
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where p 
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{Qn} and {Pk.N} are probability distributions and 



x is a two-mode pure state of a single particle [21 1 



J k,Nl 



The existence of a SSR is the consequence of the lack of 
a suitable phase reference frame (RF) [l9[. However, the 
possibility that a suitable RF can be established in princi- 
ple cannot be excluded [l9j]. If SSRs are lifted, then states 
Pcoh containing coherent superpositions of different num- 
bers of particles would become physically relevant. These 
will be called separable if they are separable in every 
fixed- TV subspace, i.e., if the state p re d = J2n 

PcohJJ-iV 
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without such coherences is of the form of Eq. ([2]) 
States which are not separable according to this defini- 
tion are entangled 22]. In this manuscript, it will be 
shown that entanglement is a necessary resource for sub 
shot-noise sensitivity in linear interferometers. 

Phase sensitivity with a linear interferometer. We con- 
sider a linear two-mode interferometer where the density 
matrix p evolves as p(9) = e~ H pe +ieJii , where 9 is 
a real number, ft is an arbitrary direction in the three 

dimensional space, and J r i — J ■ ft is a collective spin op- 
erator. For a non-fixed number of particles J is defined 

as J= ®+™J( N \ where JW = \ #0 an d |W is 
the vector of Pauli matrices acting on the Z-th particle 
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24j . Note that linear interferometers preserve the num- 
ber of particles, [ Jjj, TV] = 0. The goal is to estimate the 
parameter with the maximum possible sensitivity AO. 
The quantum Cramer-Rao theorem ensures that for ar- 
bitrary, unbiased, phase estimation protocols, the phase 
sensitivity is bounded by [25l. |26| 
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Here m is the number of independent measurements done 
with identical input states p and Fq [ J^, p] is the quantum 
Fisher information (QFI) [26j]. Equation (|3|) is a bound 
on AO given the state p and the number of measurements 
m and can be saturated in the central limit. 

In the following, we set the fundamental sensitivity 
bounds of the phase estimation problem by taking, as a 
constraint, the average total number of particles, TV to t = 
m(N), without fixing m or (TV) separately. This correctly 



10|,[27|. Aver- 



accounts for the finite available resources 
ages are computed with respect to the input state of the 
interferometer. 

Shot-noise limit. We define the shot-noise limit as 
the maximal phase sensitity attainable with separable 
states. For separable incoherent mixtures, Eq. @, 

F Q [p scp ,M = EnQnFqIp^J^} < J2nQ" n 
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(TV). The first equality is due to the block-diagonal form 
of Ja and p scp [28[ and the second inequality was proved 
in Ref. 0. By using this result and Eq. the shot- 
noise limit is thus given by 
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which agrees with the common definition of the shot- 
noise or standard quantum limit. As shown in Appendix 
A, Eq. ([4]) holds also when considering separable states 
with coherences. This brings us to the following results. 

Entanglement and sub shot-noise sensitivity. An arbi- 
trary state with non-fixed number of particles is entan- 
gled if it fulfills the inequality 



_ W 
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for some direction n. States satisfying Eq. ([5]) are useful 
in a linear interferometer implemented by the transfor- 
mation Jft, since, according to Eq. (|3]), they provide a sub 
shot-noise (SSN) phase sensitvity. We further introduce 
the spin-squeezing condition 
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< 1, 
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where the vectors fU, H2 and 713 define a right-handed 
coordinate system [29(. In Appendix B, we prove that 
X 2 < £ 2 holds. Hence, spin-squeezed states (£ < 1) are 



entangled and useful for SSN interferometry. Equations 
§5§ and © generalize the conditions for entanglement 
and SSN discussed in 0, Q to states of non-fixed num- 
ber of particles and only require the replacement of TV 
with (TV) . This result has been previously used without 
justification in several experimental works. 

The Heisenberg limit. While it is well known that 
the shot-noise limit Eq. ^ can be overcome, finding 
the Heisenberg limit (HL), i.e. the ultimate scaling of 
the phase sensitivity imposed by Quantum Mechanics 
for states with non-fixed TV, has generated a vivid de- 
bate (l. [27l. l3(il - [32j . For instance, it has been suggested 
that A9 ~ l/(A r ) is the fundamental sensitivity bound 
[3^ . Recently it was argued that the HL is given by 



A9 = 1/y (TV 2 ) 8|. Not only this limit overcomes the 

bound AO = l/(N) but, since the quantity (A^ 2 ) can 
grow arbitrarily fast with (TV), the phase sensitivity can 
be arbitrarily high when considering wildly fluctuating 
number distributions [33| . 

The definition of the HL depends on the specific con- 
straints imposed on the resources. In the presence of the 
SSR, and by taking (N), (N 2 ) and m as constraints, the 
phase sensitivity is bounded by 



AO > max 
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The first bound is derived from Eq. ([3]) by notic- 
ing that, for incoherent mixtures ((T|), FQ[p- mc , J r t} — 

ZnQnFqIpWJW] < EnQnN 2 = (N 2 ). It agrees 
with the bound of [8J, except for the factor m, account- 
ing for repeated independent measurements. Note that 
in the central limit, the bound scales as 1/^/rn. The sec- 
ond bound is proven in Appendix C. When the number 
of particles is fixed and equal to TV, Eq. (J7J recovers the 
definition of the HL AO =1/ y/mN discussed in [7| which 
takes, as constraints, TV and m, separately. Note, how- 
ever, that from this limit one does not obtain Eq. Q by 
naively replacing TV — > (TV) . 

Here, we define the HL as the maximum phase sensitiv- 
ity attainable with fixed TV to t- If coherences can neither 
be generated nor measured, as in current experimental 
intcrferometric protocols because of the SSR, the HL is 
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If the state contains coherences but only POVMs without 
number-coherences are applied, the HL is still given by 
Eq. ©, cf. Appendix C. In particular, this is true for the 
Mach-Zehnder interferometer, even if only the relative 
number of particles is measured or the parity at one exit 
port. The latter measurement was considered in Ref. [3lj |. 
By adapting the proof in Appendix C it can also be shown 
that the HL Eq. © holds for incoherent mixtures even 
if POVMs with coherences are available. 
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If the state contains coherences and POVMs with co- 
herences are used, then Eq. ((5J) is valid in the central 
limit only [34j]. Outside the central limit the bound 



A0 > 



1 
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holds, as shown in Appendix C. However, there are ar- 
guments suggesting that the HL may be given by Eq. @ 
also in this case [23, HH • 

Examples. We consider a Mach-Zehnder interferome- 
ter, corresponding to a collective rotation around the y 
axis, i.e., n = y, with a coherent state \a) a in the in- 
put mode a and a squeezed state in the input b as 
suggested by Caves [l], [l(J S|. Here a = \a\e l ^ a and 
£ = re t6 ^. When fixing the relative phase 2(j> a = 6, 
Eq. (J5J gives 



sinh 



sinh 



(10) 



For r = we have a coherent state entering one of the 
two ports and vacuum the other port. We obtain x 2 = 
1 and hence shot noise, as expected. In contrast, we 
obtain x 2 < 1 when r > and \a\ 2 7^ 0. The state is 
entangled and, for relatively small values of r, it is also 
spin squeezed, according to the definition Eq. ©. A 
direct calculation for 713 = 2; gives 

^ 2 _ (H 2 + sinh 2 r)(\a\ 2 e~ 2r + sinh 2 r) 



(|a| 2 — sinh r) 2 



In the limit \a\ 2 ^> sinh 2 r we have that £ 2 
1 for r > (and |a| 7^ 0). Spin squeezing is lost 
for \a\ 2 ~ sinh 2 r, even though the state remains en- 
tangled. Eqs (llOllip hold even for the state p rc d = 
V^r £}{a, C|1jv obtained in the presence of the SSR 

13, [37 1 . We finally notice that the state \a) Cg) \Q reaches 
the HL A6> ~ l/A tot when \a\ 2 rj sinh 2 r [ljj. 

Conclusions. In this manuscript we have extended the 
quantum phase estimation theory to the case of states 
with fluctuating number of particles. In particular, we 
have discussed entanglement and spin-squeezing and pro- 
vided the class of entangled states useful for sub shot- 
noise sensitivity. The Heisenberg limit has been de- 
fined by taking into account the finite available resources. 
Our results are relevant for most current experiments 
on the creation, detection and application of entangled 
states with a fluctuating number of particles for quantum 
metrology. 

Appendix A. Shot-noise limit. Here we demon- 
strate that Fq(p cs ) < (A) for states of the form 

Pes = 2fePfe|V'cP)(0cs I, which contain coherences be- 
tween different numbers of particles and are separa- 
ble in every fixed- A subspace. We obtain Fq(p cs ) < 

max 4X) fe Pfc(AJ f f) 2 < 4 £) fc p k Y, N max 



([4 JV) ] 2 )u,-- 



the convexity of the quantum Fisher information [38j |. 
We expanded \ipcs } 



EWQVVcs 
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for pure separable states ([Jioc' 1 ] 2 ), ,<«), < j M. By in- 
serting Fq(p cs ) < (A) into Eq. ©, we recover Eq. @. 

Appendix B. Spin- squeezing inequality. We consider, 
without loss of generality, a coordinate system such that 
( Jn 2 ) = 0. From the inequalities Fq[p, Jn] > F^^[p{0)] 

M^F m [p(9)} > (KI)^(^) 2 (obtained by extend- 
ing the proof in [391 ] to the case of non-fixed number of 
particles), it follows that X 2 < (N) (Ai) 2 /(d(A)/d<9) 2 . 
By choosing n = r?2 , A = and using the commutation 
relations of the Jft i operators [29j], we obtain iL{Jft 3 ) = 



and used that 

N 



iTr[J n3 [j H2 ,p(6)]] = (J Hl ). Hence X 2 < 



Appendix C. Heisenberg limit. In the presence 
of the SSR, only POVMs without coherences are 
available. A POVM is of this form if we can write 
its elements as E{X) = E(w M)eI(X) En,Mj where 
Em^.M = ^-N {M is a degree of freedom in each 
A-subspace), and I(X) are all pairs (A, M) leading to 
the same A (A could be the number of particles at one 
port of a Mach-Zehnder interferometer, for instance). 
For m — 1, the conditional probabilities are P(X) — 
E iN ,M )m x)P(N,M\6) = E { N,M)ei(X)QNP{M\6,N), 

where P(M\6,N) = Tr [E NjM e 



i9J H £(N) e -i8J s -\ 



and QnP^ — ^-nP^-n- The variance of an es- 
timator # C st(A) (assumed to be unbiased, i.e., 



is (A0 e 



£ A P(A|0)(0 OBt (A) 



ExE( N ,u)ei(X)QNP(M\e,N)(e est (X) - 6) 2 

lZ N Q N E M P( M \ N ^)( d e S t(x[N,M}) - ef 

EnQn^O^) 2 > J2nQn/N 2 . We used that 
E A E(iv,M)e/(A) = EatEa/ holds since Ea^( a ) = !■ 
In the last inequality we used that (A0g s1 /) 2 > 1/N 2 
holds for unbiased estimators Q. For m > I, 
P(X\6) = Y\T=iP{X l \e) and it can be shown that 

(A0 es t) 2 > E^/E™i^ > E^Qa?/(E"i^) 2 , 
where the first inequality follows as above and the second 
inequality holds for positive numbers. For m — 1, we 
then obtain (£ w Qat/A 2 )(£ w Q n ) > (E n Qn/N) 2 
from the Cauchy-Schwartz inequality. Applying it again 
leads to (EnQn/N)(J: n Q n N) > (EnQn) 2 = 1. 
Therefore, we arrive at A6* cst > 1/(A). This can be done 
in analogy for m > 1. This proves that A9 cst > l/m(N). 

The bound ([9]) can be proven as in Appendix A by 
using general states and Fq[\^ n >), J^] < A 2 . 

Acknowledgements. We thank J.I. Cirac, K. M0lmcr, 
A.S. S0rensen and T. Rudolph for stimulating discus- 



4 



[1] CM. Caves, Phys. Rev. D 23, 1693 (1981). 

[2] B. Yurke, S.L. McCall and J.R. Klauder, Phys. Rev. 
A 33, 4033 (1986); M. Xiao, L.A. Wu and H.J. Kim- 
ble, Phys. Rev. Lett. 59, 278 (1987); P. Grangier, R.E. 
Slusher, B. Yurke and A. LaPorta, Phys. Rev. Lett. 59, 
2153 (1987); M.J. Holland and K. Burnett, Phys. Rev. 
Lett. 71, 1355 (1993). 

[3] D.J. Wineland, J.J. Bollinger, W.M. Itano, and D.J. 
Heinzen, Phys. Rev. A 50, 67 (1994). 

[4] M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993). 

[5] A. S0rensen, et al, Nature 409, 63 (2001). 

[6] L. Pezze and A. Smerzi, Phys. Rev. Lett. 102, 100401 

(2009) . 

[7] V. Giovannetti, S. Lloyd and L. Maccone, Phys. Rev. 

Lett. 96, 010401 (2006). 
[8] H.F. Hofmann, Phys. Rev. A 79, 033822 (2009). 
[9] A. Rivas and A. Luis, Phys. Rev. Lett. 105, 010403 

(2010) . 

[10] L. Pezze and A. Smerzi, Phys. Rev. Lett. 100, 073601 
(2008). 

[11] Y.P. Huang and M.G. Moore, Phys. Rev. Lett. 100, 
250406 (2008); U. Dorner, et al, Phys. Rev. Lett. 102, 
040403 (2009); D.W. Berry, et al, Phys. Rev. A 80, 
052114 (2009); 

[12] P. Walther, et al, Nature 429, 158 (2004); M.W. 

Mitchell, et al., Nature 429, 161 (2004); T. Nagata, et 

al., Science 316, 726 (2007). 
[13] V. Meyer et al., Phys. Rev. Lett. 86, 5870 (2001); D. 

Leibfried et al, Science 304, 1476 (2004). 
[14] J. Appel, et al., PNAS 106, 10960 (2009); M.H. Schleier- 

Smith, I.D. Leroux, V. Vuletic, Phys. Rev. Lett. 104, 

073604 (2010). 

[15] G.B. Jo et al, Phys. Rev. Lett. 98, 030407 (2007); J. 
Esteve, et al, Nature 455, 1216 (2008). 

[16] P. Bohi, et al, Nat. Phys 5, 592 (2009). 

[17] G.C. Wick, et al, Phys. Rev. 88, 101 (1952). 

[18] B.C. Sanders, S.D. Bartlett, T. Rudolph and P.L. Knight, 
Phys. Rev. A 68, 042329 (2003); K. M0lmer, Phys. Rev. 
A 55, 3195 (1997); 

[19] S.D. Bartlett et al. Rev. Mod. Phys. 79, 555 (2007). 

[20] R.F. Werner, Phys. Rev. A 40, 4277 (1989). 

[21] Density matrices p( ' of N particles with 2 degrees of 
freedom act on a Hilbert space Hn = ®f =i H%\ where 
W]p ~ C 2 . The operator Ijv is the projector onto tLn ■ 

[22] The definition of separability can also be applied for- 
mally when treating indistinguishable particles in first 
quantization, where the set of admissible quantum states 
is reduced to the symmetric/antisymmetric subspace for 
bosons/fermions. Separable states of bosons without co- 
herences are then given by Eq. © with \cj>f N ) = \4> kjN ) 
for all i. It has been argued that particle entanglement 
due to (anti-)symmetrization in first quantization is un- 
physical since the individual particles cannot be ad- 
dressed [23| . However, in interferometers using only col- 
lective operations, it can still be a useful resource for sub 
shot-noise sensitivity. 

[23] P. Zanardi, Phys. Rev. A 65, 042101 (2002); A. P. Hines, 
R. H. McKenzie, and G. J. Milburn, Phys. Rev. A 67, 
013609 (2003); F. Benatti, R. Floreanini, and U. Mar- 
zolino, Ann. Phys. 325, 924 (2010). See also L. Amico 



et al, Rev. Mod. Phys. 80, 517 (2008) and References 
within. 

[24] For bosons, the operators J x , y ,z can be written as J x = 
(a t S + & t a)/2, j y = (a f 6 - b t a)/2i, and j z = (a) a - 
b'b)/2 using annihilation operators a and b for the two 
interferometric modes. 

[25] S.L. Braunstein and CM. Caves, Phys. Rev. Lett. 72, 
3439 (1994). 

[26] To estimate the parameter 9, a positive operator value 
measurement (POVM), {E(e)}, is performed on the 
output state and 9 is inferred from the result of this 
measurement. The maximum sensitivity of this estima- 
tion is bounded by the Cramer-Rao bound A#cr = 
1 / JrnFfe^ \p(P)\, where m is the number of measure- 
ments, F 6(e) [p(9)] = J deP(e\6) (d log P(e\9)/de) 2 is the 

Fisher information and P{e\9) = Ti[E(e)p(9)]. The QFI 
is the maximum Fq[p(9)] = max^ (en Fg u) [p(9)] and it 
is saturated by a particular POVM [25| . 

[27] S.L. Braunstein, A.S. Lane and CM. Caves, Phys. Rev. 
Lett. 69, 2153 (1992); A.S. Lane, S.L. Braunstein and 
CM. Caves, Phys. Rev. A 47, 1667 (1993). 

[28] The QFI for a state of fixed number of particles, 

is F Q [p {N \ f a N) ] = 4(A^ (JV) ) 2 , where the Hermitean 
operator R (N) satisfies the equation {R m ,p (N) } = 
i[Ji N \p W ]. For pure states, F Q U),j H ] = 4(AJ H ) 2 ^ ) . If 
the input state is an incoherent mixture p = ~^2 N QnP , 
then Fq[p, J a ] = J2 N QnFq[p (n) , jl N >] follows since p 
and Jn commute with N, which can also be imposed 
for R. The POVM saturating F = Fq can be cho- 
sen as a von Neumann measurement on the eigenstates 
of R [25|], therefore a POVM without coherences satu- 
rates Fq . However, for non-restricted states and POVMs 
Fq(p) > Fq(^2 n Ijvpljv) holds. A simple example where 
the inequality is strict is (\y+) + \y—} ® \ z +))/V2 for the 
generator J y , where \(y, z)±) are the eigenstates of <J y , z . 

[29] For a right-handed coordinate system fti, i = 1, 2, 3, the 
operators J f t i satisfy the SU(2) algebra of angular mo- 
mentum: [J^i, Jnj\ = iei,j,kJti k , where £i,j,fc is the Levi- 
Civita tensor. 

[30] Z.Y. Ou, Phys. Rev. A 55, 2598 (1997). 

[31] P.M. Anisimov, et al, Phys. Rev. Lett. 104, 103602 
(2010). 

[32] J.H. Shapiro, S.R. Shepard and N.C Wong, Phys. Rev. 
Lett. 62, 2377 (1989). 

[33] This happens, for instance, for the state proposed by 
Shapiro, Shepard and Wong (SSW) [2^]. However, the 
SSW result did not take into account the number of mea- 
surements. Indeed, it was shown that the SSW proposal 
does not overcome the limit A9 = l/N to t [27j . 

[34] Fisher's theorem ensures that the bound l/(m(7V r2 )) 1/ ^ 2 
can be saturated in the central limit for m > m c i, where 
m c i depends on the POVM used. In particular, it can be 
saturated with POVMs without coherences. In this case, 
and for m > mj,f, l/(m(A> 2 }) 1/2 > l/m(N) holds due 
to Eq. (0. However, for m' < mj." c , we cannot rule out 
that A9 < 1/iVtot can be achieved if states and POVMs 
with coherences are available. 

[35] L.B. Levitin and T. Toffoli, Phys. Rev. Lett. 103, 160502 
(2009). 

[36] The state \a) (g) |£) ((g) is the mode tensor product here) 



5 



reaches a SSN sensitivity A8 = ej \J\ct\ 2 , with e < 1 p[. 
[37] For the considered states written as \ip) = 

J2n VQiv|''/ ,(JV) )i one can show that {Ji% )i^i.n)\ = 
holds for all N. Eqs (|10lll|) are invariant since Fq(\iP}) = 
4(AJ H )^ ) = F Q (p rcd ) = E JV Qiv(A4 JV) )^(«)) H 

holds and for block-diagonal operators (-A)|^> = {A)p Ted . 
[38] Convexity of F has been proven in M. L. Cohen, IEEE 
Transactions on Information Theory 14, 591 (1968). 
A slightly more compact proof is reported here. For 



P(0) = T, k Pkpk(6) we have P(e\0) = Tr[E(e)p(0)] = 
J2 k PkPk(e\e) where P k (e\6) = Tr[E(e)p k (9)]- Us- 
ing the Cauchy-Schwarz inequality we have that 
(dP(e\6)/d8) 2 /P(e\e) < J2 k p k (dP k (e\e)/d6) 2 /P k (e\8). 
Integrating over de we obtain that F^,,-.\^2 k pi,pk(d)] < 
^2 k PkFj S ^- ) [p k (9)]. This inequality holds for all possible 
POVM and, in particular, for those saturating the QFI. 
[39] H. Uys and P. Meystre, Phys. Rev. A 76, 013804 (2007). 



